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For Anderson tight-binding models in dimension d with random on-site energies e^? and critical 
long-ranged hoppings decaying typically as V*^''{r) ~ V/r''', we show that the strong multifractality 
regime corresponding to small V can be studied via the standard perturbation theory for eigenvectors 
in quantum mechanics. The Inverse Participation Ratios Yq{L), which are the order parameters of 
Anderson transitions, can be written in terms of weighted Levy sums of broadly distributed variables 
(as a consequence of the presence of on-site random energies in the denominators of the perturbation 
theory). We compute at leading order the typical and disorder-averaged multifractal spectra Ttyp{q) 
and Taviq) as a function of q. For q < 1/2, we obtain the non-vanishing limiting spectrum Ttyp{q) = 
Tav{q) ~ d{2q — 1) as V ^ 0^ . For q > 1/2, this method yields the same disorder-averaged 
spectrum Tav(q) of order 0{V) as obtained previously via the Levitov renormalization method by 
Mirlin and Evers [Phys. Rev. B 62, 7920 (2000)]. In addition, it allows to compute explicitly the 
typical spectrum, also of order OiV), but with a different (/-dependence Ttyp{q) ^ Tav{q) for all 
q > qc = 1/2. As a consequence, we find that the corresponding singularity spectra ftyp{a) and 
fav{oi) differ even in the positive region / > 0, and vanish at different values a^^^ > a+"j in contrast 
to the standard picture. We also obtain that the saddle value atyp{q) of the Legendre transform 
reaches the termination point a^_^'' where ftyp{a^^) = only in the limit q — >■ -|-oo. 

I. INTRODUCTION 

Since its discovery fifty years ago [l|, Anderson localization has remained a very active field of research (see 
the reviews The order parameters of Anderson transitions are the inverse participation ratios (I.P.R.) of 

eigenfunctions V'(^) on ^ finite volume L'^ 



/^,d^r1V.(r)p ^ 



(the denominator can be omitted if the eigenfunction has been normalized with d'^'r\^{f)\'^ = !)• As L — > +oo, 
these I.P.R. converge to finite values in the localized phase, and behave as L'^^^^'?) in the delocalized phase. At 
criticality, the eigenfunctions become multifractal and the I.P.R. involve non-trivial exponents 0, It is actually 
useful to introduce both the typical and the averaged exponents W\ 



' L— >oo 

YjL) ^ (2) 
The distribution Pq{yq) of the rescaled variable yq — Yq{L)/Y*yP{L) is expected to decay as a power-law [9l.[l0j 



so that the typical and averaged exponents coincide or not according to the value of (3q [8|, |9| 

Ttypiq) = Taviq) if > 1 

nypiq) + Taviq) if /S? < 1 (4) 

In particular in the region q > 0, one expects that there exists a critical value Oc where /3g^ = 1 separating the region 
where the two coincide for q < qc from the region where the two differ with [1,1^ 

Pq>q. = J (5) 

For the usual short-ranged Anderson tight-binding model in finite dimension d, one expects a continuous inter- 
polation @, [ll| between a 'weak multifractality' regime obtained in the d = 2 + e expansion [T3| (the leading order 
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corresponds to the Gaussian parabolic approximation of the multifractal spectrum) and a 'strong multifractahty' 
(SM) that occur in high dimension, with the following limiting form as c? — >■ +00 

rf,f (g) = r!^^\q) = d{2q -I) iov q<\ 

4fiq)^r!^^'iq) = for g > i (6) 
For Anderson tight-binding models with long-ranged hoppings with typical asymptotic decay 

V*yP{r) oc — (7) 

r— )-+oo t'^ 

one has also found a continuous family of the critical points at a = d as a function of the amplitude V [Ijl, that 
interpolates between 'weak multifractahty' for large V — > +00 and 'strong multifractahty' for small y — >■ Q. In 
dimension c? = 1, corresponding to the ensemble of L x L power-law random banded matrices (PRBM), various 
properties of these Anderson transitions have been studied in detail [g. Il9l-^9| . In particular, the regime V ^ 
has been analyzed via the powerful Levitov renormalization method [3, |30|, that allows to compute exactly the 
disorder-averaged spectrum in the region q > 1/2 

where b = V/W is the ratio between the amplitude V and the width W of the on-site random energies. The limiting 
'strong multifractahty' spectrum of Eq. [6] has been then obtained in by using the symmetry of the multifractal 
spectrum (T3l - [l7| connecting exponents with g < 1/2 to exponents with q > 1/2 

Tav{q) - Tav{l - q) = d{2q - I) (9) 

which is expected to hold for any Anderson transition in the so-called 'conventional symmetry classes' Q. The 
derivation of Eq. |6]is thus rather indirect since one performs an explicit calculation in the region q > 1/2 to obtain 
that the exponents vanish Tav{q > 1/2) — > in this region, and the non- vanishing results in the region q < 1/2 are 
then entirely based on the symmetry of Eq. [HI The same methodology has been followed for related hierarchical 
models [3l[ and for the analysis of short-ranged models in high d [Tlj . 

The aim of this paper is to reconsider the regime of 'strong multifractahty' V ^ for Anderson tight-binding models 
with long-ranged hoppings in order to compute explicitly both the typical and the disorder-averaged multifractal 
spectra in both regimes < g < 1/2 and q > 1/2. The paper is organized as follows. In section |IT1 we describe the 
model and apply the standard perturbation theory of quantum mechanics to obtain the leading order expression of 
the I.P.R. Yq of Eq. [T] In section IIIIl we analyze the statistical properties of some weighed Levy sums that play a 
major role in this perturbation theory. The typical behaviors and the disorder-averaged behaviors of the I.P.R. Yq are 
then computed in sections IIVI and |V] respectively. Our final results concerning the typical and averaged multifractal 
spectra in the various regions of q are discussed in section IVTl and are compared to previous results. Our conclusions 
are summarized in section fVIII 

II. PERTURBATION THEORY FOR ANDERSON TIGHT-BINDING MODELS WITH 

LONG-RANGED HOPPINGS 

A. Anderson tight-binding models with long-ranged hoppings 

We consider an Anderson tight-binding model on an hypercubic lattice of size L'' defined by the Hamiltonian 

Ho ^ 

r 

Hi = J2 yp,P'\r><r'\ (10) 

r^r ' 

where the hoppings are symmetric / = V^*^^ / with = 0. 
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1. Assumptions on the random on-site energies ep 

We consider the usual case where the on-site energies are independent random variables distributed with a law 
P(e) which is symmetric around e = and which presents the scaling form 

n^) = ^rp{^) (11) 



so that W represents the disorder strength. In the following, an important role will be played by the probability 
density around e = that will be denoted by 

P(. = 0) = ^ (12) 

where c = 2p (0) > 0. For instance in the PRBM model [8,], P(e) is a Gaussian of variance unity corresponding to the 
values 

^PRBM ^ ^ 

c™^^^ = (13) 
The average over the random energies {e^} will be denoted by 



E{Ai{e,})) 



A{{er}) (14) 



2. Assumptions on the long-ranged hoppmgs Vp^ft / 
We consider the critical case where the long-ranged hoppings decay typically as in dimension d [ll, [l^ 

V{r} ^ (15) 

r— 7-00 r 

where up are either fixed {up = 1) or are independent random variable of order 0(1). For instance in the PRBM 
model in c? = 1 Q, itp is a Gaussian of variance unity. In any case, we assume here that the distribution of up is such 
that its moments exist. 

In the following, two quantities will play a major role. Denoting by Sd the surface appearing in the radial change 
of variables d'^r = Sdr'''~^dr in dimension d, we may evaluate the behavior in L of the following sums 

Y,WW)\ = SdVW\\nL (16) 

r 

and for q <l/2 

2n f 1 rd{l-2q) 

Y: \Vir-)\ = V'^\up\^^ J Sdr'-'dr-^ = v'^\up\^^Sd—^ (17) 

r 

B. Perturbation theory in the hoppings 

As explained in the Introduction, we focus in this paper on the 'strong multifractality regime' that corresponds to 
a small amplitude V in the long-ranged hoppings of Eq. 1151 It is thus natural to consider the perturbation theory 
associated to the decomposition of Eq. [101 The Hamiltonian Hq has for cigcnstates the L'^ completely localized 
eigenfunctions on each lattice site, the eigenvalues being simply the corresponding random on-site energies ep 

l'/'f> =\r> 

4°' = (18) 



4 



The standard perturbation theory of quantum mechanics yields that at lowest order, the eigenvalues are unchanged 

= (19) 

whereas the eigenfunctions read 



W > 

The corresponding I.P.R. of Eq. [T]read 



i# > + E > 



a(0) 



> 



r > 



r ' > 



(1) _ 



Y.r'\€\r')\ 



in terms of the sums 



(l + Si)'' 



29 



(20) 



(21) 



(22) 



To simplify the notations from now on, we will focus on the eigenstate associated to the central site r — Q and we will 
consider that its associated eigenvalue is exactly at the center of the band eg = 0. Then the perturbed eigenenergy of 

Eq. [19] is also i?!^"* = ejj = at leading order, and the I.P.R. of the corresponding eigenstate (Eq. [2T|) will then allows 



to compute at leading order the multifractal spectrum corresponding to = 0. The variables of Eq. 



V{r) 



fir 



2q 



] become 
(23) 



The aim of this paper is to analyze the statistical properties of the I.P.R. of Eq. [5TJ in particular their typical values 
and their disorder- averaged values to extract the multifractal exponents Ttyp{q) and Tav{q) of Eq. [2l It is convenient 
to analyze first the statistics of the sums Eg that turned out to be weighted Levy sums as we now explain. 



III. 



STATISTICS OF THE WEIGHTED LEVY SUMS S„ 



In this section, we discuss the statistical properties of the sums of Eq. 1231 that can be rewritten as sums 

S,^ ^ \Vif)\"^z,if) 



of the random variables 



(24) 



(25) 



with the weights \V{f)\'^''. 

We recall here that the average of an observable O over the random on-site energies will be denoted by E{0) (see 
Eq. I14p . In the case where the long-ranged hoppings of Eq. [T5]are non-random (up = 1), the disorder-average denoted 
by O is equal to E{0). In the case where the long-ranged hoppings of Eq. [15] are also random, the disorder-average 
denoted by O denotes the average over both the random on-site energies and the random variables up appearing in 
the long-ranged hoppings of Eq. [T51 



Statistics of the variables Zqif) = |e,^ 



-2q 



From the probability density P{e ~ 0) near zero energy given in Eq. 1121 one obtains via a change of variable that 
the probability distribution Qq{zq) of the variable Zq(r) = \ep\~'^'^ presents the following power-law decay 



Qqi.Zq) ^ 



with n„ = — 
2q 



(26) 
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In particular, the disorder-averaged value E{zq) (with the notation of Eq. [14)) presents a transition at g = 1/2 

1 

2 



E{zq) < +00 for q > - 



E{zq) = +0O for q<- (27) 



B. Generating function E{e of the sum 

Levy sums of identically broadly distributed variables (without the weights |l^(r)p'' in Eq. ^IM appears in various 
fields of disordered systems, usually in low-temperature disorder-dominated phases of classical models : their statistical 
properties are described in particular in [36 39]. In the following, we analyze the effects of the presence of the weights 

|y(f)p^. 

The first important property is that the disorder-averaged value E (S^) presents the same phase transition as Eq. 
[27l as a consequence of the linearity of the sum of Eq. [24l 

1. Case q < 1/2 where the disorder-averaged value E {T,q) is finite 

For q < 1/2, the disorder- averaged value is finite and reads 

Eii:,) = Eiz,) ^ iVir^'^ (28) 

where using Eq. [TT] 

/+00 r+oo 
deP{e)\e\''^'^ ^W-^Wg with Bg ^2 dxp{x)x-^'^ (29) 
-oo Jo 

For instance if p{x) is Gaussian of variance unity as in the PRBM model, one obtains 

/•+°° p~4- Vi- — a\ 

^Gauss ^ 2 / da;^=a;-^« = ^' (30) 
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After averaging also over the possibly random hoppings of Eq. 1151 we finally obtain using Eq. [17] 



2. Case q> 1/2 where the disorder-averaged value i5(Eq) is infinite 

For g > 1/2 where < 1, the generating function E (e^^^i) presents the characteristic singularity in t^i of Levy 
distribution (see Eq. [26]) 

E ^ dzMzq)e-'^-^ = ^ ' dzqQq{zq){l - e"*^') 1 - ^ [-r(-M,)] + - 

e-*"'^[-r(-M,)]+... (32) 



t-fO 

with the usual integral 







The generating function E (e"*^?) will thus presents a similar singularity 



+ 00 7 

T^(l-e--) = -r(-^,) (33) 



~ e-*'''^["^("^''l^''l^(*''l (34) 
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After averaging also over the possibly random hoppings of Eq. [151 '^e finally obtain using Eq. [TBI 



~ g-t^''^[-r(-/.<,)]|«,HS<ilnL /ggN 
t->-0 



Equivalently, inverting Eq. 1321 one obtains the following asymptotic decay for the probability distribution Vq(T,q) 



VJT,,) ~ — r-S- (36) 



with the exponent /ig = l/{2q) and the amplitude 

F 

Aq = cSdfiq — \uf\\nL (37) 

C. Analysis of the auxiliary quantity E (^'e^'^^) 

The variables and Si that appear in the numerator and denominator of Eq. [5T]are correlated since they involve 
the same random energies. In the following, some computations will involve the auxiliary quantity 

E{zle-''') (38) 

From Eq. 1261 the probability density Qi{zi) is know to decay with a power-law of exponent 1 + jJi ~ 3/2 

Qi(^i) ^3 72 (39) 

As a consequence for q < 1/2, the non-integer moment E{z1) of order q exists and one has 

E{zte-'^^)^^^Eiz'[) + o{t) (40) 

For q > 1/2, the non-integer moment E{z'[) of order q does not exist, and the divergence of the auxiliary quantity as 
t — >■ can be obtained via the change of variables x = tzi 

/■+°° clr /x\ /x\<i 



E{zle-^^-) ^ ^ dzMz,)zle-'^^ = °° ) - e 



C 1 „ /■+°° , 3 „ _^ C 1 / 1 



-^5-9 / dx^-'^xe-'' = r^-iT [q--] (41) 



To obtain the typical values of the I.P.R. Yq (see Eq. [5]), we need to average the logarithm of the expression of Eq. 



lni;*^P = \nYq = In (1 + E,) - gin (1 + Si) (42) 

So here the correlations between Sq and Ei do not play any role, and one only needs to know the statistical properties 
of the sums and Ei discussed in the previous section. 



A. Computation of ln(l -I- E,) for g < 1/2 

For q < 1/2, the disorder-averaged value Eg converges (see Eq. [3T|) . So we may use the expansion of the logarithm 
to obtain the leading-order behavior 

Ld(l-2q) 



ln(l + E,) BqS,W\-^^^^^ (43) 



B. Computation of ln(l + Eg) for q > 1/2 

For q > 1/2, we use the following integral representation of the logarithm 

ln(l + S,) = ye-*(l 



to relate the disorder-averaged value to the generating function e *^<j 



ln(l + !],)= / ^e-* (l - e 



Using the result of Eq. [35l we obtain 



ln(l + Sq) 



Using the relation 



the final result reads 



r(-Mg)r(Mg) 



sin(7r^q) 



ln(l + E,) ~ Si\x,L—\uf\- 



(44) 



(45) 



(46) 
(47) 

(48) 



W sin(7r/Xg) 

Note that this result can be equivalently obtained via a direct calculation based on the asymptotic behavior of Eq. 



ln(l + E,) dE,7',(S],)ln(l + E,)~^, / ln(l + S,) = X 

-'0 -'0 Eo 



/i, sin(7r/i. 



= 5'rflnLlufl 



CU TT 

sin(7r^q) 



(49) 



C. Behavior of the typical I.P.R. V, for < g < 1/2 



For q < 1/2, we have seen that ln(l + E^) grows as L'^^^ fcf |43|). whereas ln(l + Ei) grows only as InL (Eq 
. We thus obtain at leading order 



W ) 



Ld{l-2q) 

d{l - 2q) 



The typical exponents Ttyp{q) defined in Eq. [2] thus read 

Ttyp{q< 1/2) ^d{2q-l) 
in agreement with the 'strong multifractality' limit of Eq. [6l 



(50) 
(51) 



D. Behavior of the typical I.P.R. Y, for g > 1/2 



For q > 1/2, we have seen that both ln(l + Eg) and ln(l + Ei) grow as InL (Eq. |48|) . These two contributions 
yield at leading order 



cUt 



InFg =ln(l + Eg)~gln(l + Ei) = Sd\nL—\up\ 



sin(7r/ig) '^sin(7r/zi) 



Since /xi — 1/2, we obtain that the typical exponents Ttyp{q) defined in Eq. [5] read 



cV 

Ttyp{q > 1/2) = Sd—\uf^\Tr 



sin(^) 



(52) 



(53) 



V. DISORDER-AVERAGED VALUES OF THE I.RR. V, 

To compute the disorder-averaged values of the I.P.R. Yq of Eq. [211 it is convenient to use the identity 

1 1 



a« r(g) Jo 

to obtain 



+ 00 



— -^l/irsf: contribution H~ second contribution 

1 /' + ^ _^ 

first contribution — TT7 7 / 6 C ^ ^ 

r(g) 7o 

T /• + °° 

^ q\secona contribution -r^/ \ / ^q"^ 

r(g) 7o 

We now evaluate separately these two contributions. 

A. Computation of the first contribution 
Using Eq. [35] for g = 1 with /ii = 1/2 and [-r(-l/2)] = 2^ 



the first contribution of Eq. [55] reads at leading order 

I o 



nq)Jo nq)Jo V 

1 , ^ „, Mr lgrflnL 



r(«) 

B. Computation of the second contribution 

To evaluate the second contribution of Eq. [55] we first need to evaluate, using the definitions of Eas[24land[ 

r 

For g < 1/2, we use Eqs [40] and [56] to obtain 

r 

SO that the second contribution reads at leading order using Eq. [IT] 

J^d(l-2q) 

Yg<if2\secondcontribution ^ ^ E (z?) ^ \V {f)^''' ~ E (zf) ^rfT/'^ I^H''^ _ ^ n 

r 

For g > 1/2, we use Eos [4T] and [56] to evaluate the singularity of Eq. [58] for small t 
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Using Eq. 1161 we finally obtain at leading order 

C. Disorder-averaged I.P.R. for q < 1/2 

For q < 1/2, the first contribution of order In L (Eq. [57| is negligible with respect to the second contribution of Eq. 
[6OI which leads to 

Ld(l-2q) 



The disorder-averaged exponents Tav{q) defined in Eq. [2] thus read 

Taviq < 1/2) = d(2g - 1) = Ttyp{q < 1/2) (64) 
and coincide the the typical exponents Ttyp{q), in agreement with the 'strong multifractality' limit of Eq. [SI 

D. Disorder-averaged I.P.R. for q> 1/2 

For q > 1/2, we add the two contributions of order InL obtained in Egs [57l and |62] Using r(z + 1) = zT{z), this 
leads to 



Cy^SdV- r 

\up\ in L 

2W ' ' 



= ^' r(g-i) l"-^'^""- (^^^ 

The disorder-averaged exponents Tav (<z) defined in Eq. [5] can be thus identified by the following expansion 

Yq(L) ^ L-^^M = 1 _ ^^^(^) (66) 

This yields 

r (o > i) = ^^ii^^^^^M (67) 

2^ r(g-l) W ' ^ ^ 

This expression coincides with Eq. [8] obtained previously via Levitov renormalization [9|] for the PRBM model with 
the correspondence b — V/W if one uses c — (Eq. [T3]) and 

e--|u| = J (68) 



-oo V27r V TT 

with d = 1 and Si = 2 (to take into account the ring geometry). 

VI. SUMMARY OF THE RESULTS AND DISCUSSION 

In this section, we summarize and discuss the results obtained for the typical and the averaged multifractal spectra 
starting from the perturbation formula of Eq. 1211 



10 



A. Typical and averaged multifractal spectra Ttyp(q) and Tav[q) 

In the region q < 1/2, we have found that the typical and the averaged multifractal spectra coincide (Eos [STl and 
I64p and are given by the 'strong multifractality' limit of Eq. [5] 

Ttypiq < 1/2) = Tay{q < 1/2) = d{2q - 1) (69) 

In the region g > 1/2, we have found that the typical and the averaged multifractal spectra are different (Egs [53l 
and 1171) 



cV 

Ttypiq > 1/2) ^ Sd—\uf\Tr 



1 



TaAq> 1/2) ^ S,^\up\V^y^-^ (70) 

where the results for Taviq > 1/2) coincides with the result of Eq. |S] obtained previously via Levitov renormalization 
i- 

Our conclusion is thus that the critical value qc where the typical and averaged spectra separate (Eq. |4]) is 

qc = l (71) 

in contrast with the other value qc ~ 2.4 predicted in [9], on the basis of the vanishing of the averaged singularity 
spectrum /at, (a) (see below Eq. \ET\i . Moreover, we expect that the distribution Pq{yq) of the rescaled variable 
Hq = Yq{L) /Y^yP{L) will decay in the scaling regime with the power law of Eq. |3] of exponent (1 + f3q) that will 
coincide with the exponent (1 + /iq) describing the distribution of T,q (see Eq. I36p 

in agreement with Eq. [5l 

B. Consequences for the typical and averaged singularity spectra ftyp{ct) and fav(ct) 

Besides the multifractal exponents Ttyp{q) and Tav{q), it is usual to introduce the typical and averaged singularity 
spectra ftyp{o) and favi^) which describe the numbers A/'^^^''^^ (a) of points r in a sample of size L'', where the weight 
|V'L(r)p scales as L"" [| 

A/'f (a) oc L-^*— (73) 
The I.P.R. Yq can be then rewritten as integrals over a 

YtyP'^^(L) ~ [ da L/'fp.-(") L-?" ~ i-Ttyp.aAq) (74) 

^ J L^oo 

The exponents Ttyp av (<?) can be obtained via a saddle-point calculation in a to obtain the Legendre transform formula 
ii 

- Ttyp,av{q) = max [ftyp,av{a) - qo] (75) 

a 

At leading order, the 'strong multifractality' limit of Eq. [HK or Eq. [52] above) corresponds to the singularity spectra 

i 

/typ(a) = /a.(a) = I for < a < 2d (76) 

The typical exponent atyp where ftypicttyp) = d thus corresponds to the maximal value atyp = otmax = 2d. The 
singularity spectrum vanishes only at the other boundary amin = 0. From the point of view of the Legendre transform 
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FIG. 1: (Color on line) (a) Comparison of the typical and averaged multifractal spectra Ttyp{q > 1/2) and Taviq > 1/2) 
of Eq. [77] : they are close near q — >■ (1/2)+ (Eq. I78[) but are very different at large q (Eq. [79]). (b) Comparison of the 
corresponding typical and averaged singularity spectra Ftyp{A) and Fav{A) : Ftyp(A) exactly terminates at the point A^^'' 
of Eq. [S2] corresponding to q — +oo, where it vanishes -Ftyp(^+^) = 0, whereas FaviA) vanishes at another value A^j^ 
corresponding to q°^ (see Eq. [8T]) and continues in the negative domain FaviA) < 0. 



of Eq. [TS] this case is singular since the saddle a*{q) is concentrated on the two values a*{q < 1/2)) = a^ax = 2c? 
and a*{q > 1/2)) = q;™„ = 0. 

Let us now take into account the small corrections of Eq. [70] in the domain q > 1/2, where we factorize the small 
prefector e = Sd^lMflf (this choice of constants in e has been preferred to have the same normalization of Tav{q) as 



i) 



Ttyp{q > 1/2) 



Ta,{q > 1/2) 



Ttypiq > 1/2) 



r„,„(g>l/2) _ 2r(q-i) 



0Fr(g - 1) 



(77) 



The two multifractal spectra Ttyp{q > 1/2) and Tav{q > 1/2) are shown on Fig. [T](a) for comparison : there are very 
close in the region q — >■ 1 /2 where they present the same singularity 



Ttypiq > 1/2) 



Taviq > 1/2) 



1 



g^(i/2)+ 7r(g-i) 
1 



9^(1/2)+ 7r(g-i) 

but are completely different at large q with the following asymptotic behaviors 

2~ 



(78) 



Ttypiq > 1/2) 



2 1 



„l/2 



Taviq > 1/2) ^ 

As discussed in Q , the disorder-averaged singularity spectrum fav (a) takes the scaling form 



favia) = eFav M = - 



(79) 



(80) 



where FaviA) is the Legendre transform of Taviq) ■ its properties have been described in In particular, it vanishes 
at A"^ corresponding to with the numerical values Q 



Al" ~ 0.51 
ql" ~ 2.405 



(81) 
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Similarly, the typical singularity spectrum ftypipt) takes the scaling form 

ftyM) = ^Ftyj, [a = ^) (82) 

where Ftyp{A) is the Legendre transform of Ttyp{q). From the asymptotic linear behavior of Ttyp{q) (see Eq. I79p . one 
obtains that Ftyp(A) exactly terminates at the point 

A'^P = 2 (^1 - ~ 0.727 (83) 

corresponding to g = +oo, where it vanishes Ftyp{A^fP) = 0. The two singularity spectra Fav{A) and Ftyp{A) are 
shown in Fig. [1] (b) for comparison. 

The fact that Ftvp{A) exists only in the region where it remains positive Ftyp{A) > is a standard property of 
typical spectrum [8|. What is surprising however is that the typical and averaged singularity spectra differ even 
in the region where there are positive, whereas the standard picture is that Ftyp{A) = Fav{A)6{Fav{A) > 0) 
Equivalently, in this standard picture [8|, the typical spectrum is expected to be exactly linear Tfyp{q) ~ for 
q > q+ meaning that in the Legendre calculation of Eq. [TSJ the saddle value remains frozen at Atyp{q > q+) — A+. 
Here we have found instead that the typical spectrum Ttyp{q) is not exactly linear in q in the region where it is 
different from Taviq), and that the saddle value Atyp{q) of the Legendre transform of Eq. 1751 reaches the termination 
point A*^P only in the asymptotic regime q — ^ +oo. 



C. Discussion on the method 



The perturbative calculation of the multifractal spectrum in the strong multifractality regime can only be very 
singular since one starts from a complete localized basis to construct multifractal critical eigenvectors via perturbation. 
To face this difficulty, two strategies have been proposed : 

(i) The powerful Levitov renormalization method [sol [3l| performs iterative changes of bases to take into account 
the resonances that occur at various scales. This approach has been reformulated as some type of 'virial expansion' 
in Refs [3|i|. 

(ii) In the present paper, we have proposed instead to use the standard perturbation theory of quantum mechanics. 
It is thus simpler than (i), since we work in the initial completely localized basis. However, since the random 
perturbation terms are singular, the essential point in this approach is that the I.P.R. should be computed with Eq. 
1211 where the perturbation terms of the eigenvectors appear both in the numerator and in the denominator : this ratio 
is then regular, because any potential divergence appearing in the numerator is compensated by the corresponding 
divergence in the denominator. In the present paper, we have described in detail how the first-order expression of 
the perturbed eigenvector allows to obtain the leading order of the multifractal spectrum in various regions of q. A 
natural question is how this approach can be pursued at higher orders. We stress that one should not use the standard 
Rayleigh-Schrodinger expressions for the normalized perturbed eigenvector (since these expressions are in fact based 
on the perturbative expansion of the normalization, which is singular here as explained above) . We believe that the 
correct formulation of our approach at higher orders involve the perturbative expansion of the eigenvector 

+00 +00 

id^r > = > + E i^^> + E > (84) 

n— 1 n— 1 

in the so-called intermediate normalization defined by 

<4°Vr>=<^10?*>=l (85) 
so that all corrections are orthogonal to the zeroth order term >= \r> 

< 4°'l4"' >=< ^4"^ >= for n > 1 (86) 
The corresponding I.P.R. of Eq. [1] should be then obtained as 



- q 

2 



^ (87) 

1+ E \^T{r') 
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VII. CONCLUSIONS AND PERSPECTIVES 



In summary, we have show that that the strong multifractahty regime of Anderson tight-binding models in dimension 
d with critical long-ranged hoppings can be studied via the standard perturbation theory for eigenvectors in quantum 
mechanics. The Inverse Participation Ratios Yq{L), which are the order parameters of Anderson transitions, then 
involve weighted Levy sums of broadly distributed variables, as a consequence of the presence of on-site random 
energies in the denominators of the perturbation theory. We have computed at leading order the typical and disorder- 
averaged multifractal spectra Ttyp{q) and Tav{q) as a function of q. For q < 1/2, we have found the non-vanishing 
limiting spectrum Ttyp{q) = Tav{q) — d{2q — 1) as F — > 0+, that had been obtained previously in [§] only indirectly 
via the symmetry relation of Eq. [9l For q > 1/2, we have obtained the same result for disorder-averaged spectrum 
Tav{q) at order 0{V) as obtained previously via the Levitov renormalization method [9]. This agreement between 
these two completely different approaches is a good indication in favor of the exactness of this result. But in addition, 
our present approach allows to compute explicitly the typical spectrum (that has not been computed via Levitov 
renormalization) : it is also of order 0{y) but has a different g-dependence Ttyp{q) ^ Tav{q) for q > qc = 1/2, and 
is not exactly linear in this regime, in contrast with the standard picture @. As a consequence, we have found that 
the corresponding singularity spectra ftyp{a) and fav{oi) differ even in the positive region / > 0, in contrast with the 
standard picture where they coincide in the positive region ftyp{a) = fav{ct)0{fav{<^) > 0) Q, and that the saddle 
value Atyp{q) of the Legendre transform reaches the termination point A*_^^ where ftyp{A*^^) = only in the limit 
q — +QO. 

In conclusion, the present work based on a pedestrian perturbative explicit approach thus questioned important 
statements of the standard picture of multifractahty at Anderson transitions. We hope that it will stimulate further 
studies to better understand the differences between typical and averaged multifractal spectra. 
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